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Abstract 

We investigate the random matrix configurations for two or three interacting 
electrons in one-dimensional disordered systems. In a suitable non- interacting 
localized electron basis we obtain a sparse random matrix with very long tails 
which is different from a superimposed random band matrix usually thought 
to be valid. The number of non-zero off-diagonal matrix elements is shown to 
decay very weakly from the matrix diagonal and the non-zero matrix elements 
are distributed according to a Lorentzian around zero with also very weakly 
decaying parameters. The corresponding random matrix for three interacting 
electrons is similar but even more sparse. 
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There is a great current interest in the locahzation weakening effect due to the 
interaction of two electrons in one-dimensional (ID) disordered systems 0, ^, ^ 
B 0' B Shepelyansky |^ mapped this problem onto a class of random banded 
matrices with strongly fluctuating diagonal elements, being the eigenenergies of the 
non-interacting problem, and independent Gaussian random off-diagonal matrix el- 
ements of zero average and typical strength lying in a band of width the one 
particle localization length i^i with U the strength of the interaction. Moreover, by 
the mapping to a superimposed banded random matrix ensemble [SBRM) a fraction 
of 1^ states with a considerable enhancement ^ ~ U'^^i of the localization length 
along the center of mass coordinate is predicted due to coherent propagation of the 
electronic pair. It was also shown that the interaction has no effect for the majority 
of other states with the two particles localized in isolated spatial positions which do 
not allow overlapping. This conclusion was confirmed and extended to higher di- 
mensions via Thouless block scaling picture by Imry [0. The subsequent numerical 
studies p, H, 1^, |, 1^ H, 1^, verified the main qualitative results concerning the 
presence of Shepelyansky states, mostly by supressing single particle transport via 
efficient Green function or bag model methods which examine pair propagation [§. 
The deviations from the predicted behavior of the two-particle localization length ^ 
found were usually attributed to the oversimplified statistical assumptions concerning 
the band random matrix model of the original Shepelyansky construction. 

However, there is an ongoing debate whether coherent pair propagation actu- 
ally exists for two interacting electrons in infinite disordered systems |TT|, O], which 



began by a recent transfer matrix study where no propagation enhancement is found 
aX E = for an infinite chain ||rT|. Moreover, it was pointed out that the reduction to a 
SBRM relies on questionable assumptions regarding chaoticity of the non-interacting 
electron localized states within ^i, so that the relevant matrix model could be prob- 



ably different ||T3| , p!4| . Although the reported absence of propagation enhancement 
|[TT| can be critisized, since the transfer matrix method may not measure the actual 
pair localization length, it is correct that the fraction of the Shepelyansky states will 
eventually shrink to zero when the system size increases, although not affecting their 
physical significance. A different localization length enhancement for theses states, 
of the form ^ ~ U'^^l'^'^ with a ^/-dependent exponent 7 < 1, was also proposed on 
the basis of numerical data by a reduction mechanism to another appropriate random 



matrix model ||T^. It must be mentioned that the Shepelyansky states are expected 
to exist as long as the interaction is not too large, since it is firmly established |]T0|, |13[ 
that in the strong interaction limit f/ ^ 00 no coherent propagation enhancement is 
possible with these states decoupled from the main band with ^ ~ ^1. 

It is worthwhile to check the validity of the mapping to a SBRM which al- 
lowed most of the previous results concerning the Shepelyansky states to be derived. 
The fact that a mapping to SBRM neglects phase correlations of the one particle 
localized states was also emphasised for a few known examples in another recent 
study [Q, where it lead to a non-justified propagation enhancement. In order to 



shed light on the appropriate random matrix description for the problem of two in- 
teracting electrons in a random potential we examine explicitly the structure of the 
two and three-electron Hamiltonian matrices by a direct numerical analysis. The 
corresponding Anderson- Hubbard Hamiltonian p!5|, |T6[ can be written as 

N N 

H =Y. IZ(^Ul,<xCn,<x + H.C.) + ^rA,a(^n,a 
n=l o" n=\ CF 

N 

+ ucl 

n=l a^a' 

where ^ and Cn^a are the creation and destruction operators for the electron at site 
n and with spin a, is energy level at site n which is a random variable uniformly 



distributed in the range [—W/2, W/2] as for the Anderson model and U is the strength 
of the interaction between the electrons. 

A non-interacting disordered system of atoms with f/ = in Eq. (1) has 
A^ linearly independent one-electron localized wave functions ipi with corresponding 
eigenvalues Ei, i = 1,2,3, A^. These exponentially localized wave functions are of 
the approximate form 



Mn) ^exp 



\n — 



(2) 
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where Ui is the localization center of ipi and 6i a corresponding phase factor. The 
corresponding perturbational localization length for small W is ~ where 
E is the single particle energy, although for some special E the prefactor is different. 
In order to examine the few-body problem in the presence of disorder we have ob- 
tained numerically all the localized states ipi and arranged them so that if i < j the 
coordinate for the localization center of ipi is smaller than the corresponding center 
rij of the wave function ipj. This is a natural kind of arrangement which guarantees 
the largest overlapping of the wave functions to occur when their indices have the 
smallest difference. 

Firstly, we consider two electrons with opposite spins and use the A^^ products 
of the two one-electron wave functions 

¥^^=^,^,,z,j = l,2,..,N, (3) 

as convenient basis states for the two interacting electrons. We do not consider the 
spin configuration of the two-electron wave function as suggested in Ref. [|^]. The 
index m is also arranged in such a way so that states with the smallest difference of 
their indices have the strongest coupling. One can calculate the matrix elements of 
the ^/-dependent Hubbard interaction term of Eq. (1) in the obtained new basis set 
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via 



(4) 



n 



In the original SBRM construction the matrix element Hm,m' was shown to vanish 
unless all four relevant wave functions were overlapping. Then each ipi{n) was assumed 
completely random within by taking the approximate states of Eq. (2) with a 

random phase factor 9i{n), so that one immediately obtains the estimate for 



the typical magnitude of the off-diagonal matrix elements distributed within a band 
range. The SBRM construction ignores phase correlations which are known to give 
regular or fast oscillations within for the single particle localized wave functions. 

In Figs. 1 and 2 we show the obtained distribution of the diagonal and the 
off-diagonal matrix elements for two interacting electrons with disorder extent W and 
interaction strength U. The diagonal matrix elements are seen to obey a Gaussian 
distribution as expected. The off-diagonal matrix elements are found to be mostly 
zero but also very close to zero with a Lorentzian distribution having very long tails 
in the latter case. In order to fit the obtained distribution for the off-diagonal matrix 
elements h in the adopted basis we use the following function sum 



where the Gaussian-hke term fa accounts for the distribution of the off-diagonal 
elements which are very close to zero and the Lorentzian-like term fi for the long 
tails. The variations of the obtained four fitting parameters Oi, 02, 03 and 04 as a 
function of the distance of the matrix element locations from the main matrix diagonal 
are shown to decay extremely weakly in Figs. 3 and 4. It can immediately be seen 
that ai is much larger than which implies a very large number of near-zero matrix 
elements. The obtained parameters 02, 04 characterising the distribution are found 
very small. 




(5) 
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We have also investigated the random matrix structure for the corresponding 
three electron problem. In order to simplify the calculations we have restricted our 
consideration to three electrons where two of them have spin up and one spin down. 
For a chain of A'"-sites in the adopted three-electron basis there are N'^{N — l)/2 
linearly independent states whose wave functions can be written as 

= i^iAMk,l (6) 

We again sort these 3-particle states in such a way so that those with the smallest 
difference in their indices have the strongest possible coupling. In Fig. 5 we display 
the obtained distributions for the diagonal and the off-diagonal matrix elements. For 
the diagonal elements the distribution is more sharp when compared to the two- 
electron case, reflecting a smearing effect in the total energy fluctuations due to the 
increase in the number of particles. We find that the off-diagonal matrix elements 
obey the general principles of the two-electron problem but the matrix structure is 
comparatively even more sparse in this case. 

In summary, we have studied the random matrix configurations for two or 
three interacting electrons in a disordered chain. We made a basis set arrangement 
which consists of non-interacting particle wave function products of orbitals so that 
the off-diagonal matrix elements due to the interaction occur successively only for 
states having closely spaced indices. Although this is the most favourable basis set 
for obtaining a SBRM a very sparse random matrix structure is shown to emerge 
from our data instead, having no well-defined band region with Gaussian matrix 
elements. We think finding another basis set consistent with a band random matrix 
structure is not easy and the few-body problem can be studied via extremely sparse 
random matrix configurations. It must be pointed out that the obtained sparse matrix 
is not incompatible with states having enhanced localization length along the center 
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of mass coordinate, since they can be shown by methods which do not rely on the 
specific matrix mapping p], p]7| . 
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Figure Captions 

Fig. 1. The distribution of the matrix elements for the two-electron Hamiltonian 
with disorder W = 3, interaction strength U = 1 and chain size = 100. (a) The 
diagonal matrix elements, (b) The first off-diagonal matrix elenents (continuous line) 
and the fitting curve from Eq. (5) (broken line). 

Fig. 2. The same as Fig. 1 but for W = 6 and interaction strength U = 4. 
Fig. 3. The four fitting parameters ai, 02, 03, 04 of Eq. (5) for the distribution of the 
off-diagonal matrix elements corresponding to Fig. 1 as a function of the distance from 
the main matrix diagonal, with chain size N = 100, disorder W = 3 and interaction 
strength U — 1. 

Fig. 4. The same as Fig. 3 but for 1^ = 6 and interaction strength U — A 
corresponding to Fig. 2. 

Fig. 5. The distribution of the matrix elements for the three-electron Hamiltonian 
with disorder W — 6, interaction strength U — A and chain size N — 100. (a) 
The distribution of the diagonal matrix elements, (b) The distribution of the first 
off-diagonal matrix elenents. 
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